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< I [ Abstract 

-t— ' . 
^ , We study the low-temperature thermodynamic properties of a number of frustrated quantum antiferromag- 

^ I nets which support localized magnon states in the vicinity of the saturation field. For this purpose we use 1) a 

mapping of the low-energy degrees of freedom of spin systems onto the hard-core object lattice gases and 2) an 

exact diagonalization of finite spin systems of up to A'^ = 30 sites. The considered spin systems exhibit universal 

2:^ ' behavior which is determined by a specific hard-core object lattice gas representing the independent localized 

Q ' magnon states. We test the lattice gas description by comparing its predictions with the numerical results 

O , for low-lying energy states of finite spin systems. For all frustrated spin systems considered we find a strong 

variation of the low-temperature specific heat passing the saturation field and a maximum in the isothermal 

entropy at saturation field resulting in an enhanced magnetocaloric effect. 
> 
Cn ■ PACS codes 



75.10.Jni Quantized spin models 



■^^ ' 75.45.-l-j Macroscopic quantum phenomena in magnetic systems 

r^ . 75.30.Sg Magnetocaloric effect, magnetic cooling 

o 
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^ ; 1 Introductory remarks 

rj-^ . Quantum spin antiferromagnets on geometrically frustrated lattices have attracted much attention during last 
C , years [1-3]. A new and rapidly developing direction in this field of quantum magnetism is a study of the low- 
^ , temperature properties of such systems in the presence of an external magnetic field. The Zeeman interaction 
of spins with a magnetic field competes with the frustrating antiferromagnetic interspin interactions that may 

K^ , lead to new phenomena. Recently it has been found that a wide class of geometrically frustrated quantum spin 
antiferromagnets (including the kagome and pyrochlore antiferromagnets) has quite simple ground states in the 
vicinity of the saturation magnetic field [4-6] (for a review, see Refs. [7, 8]). These ground states consist of 
independent (i.e. isolated) localized magnons in a ferromagnetic environment. The localized magnon states were 
used to predict a ground-state magnetization jump at the saturation field [4-6], a magnetic field induced spin-Peierls 
instability [9,10], and a residual ground-state entropy at the saturation field [3,8,11-13]. Moreover, in Refs. [8,12,13] 
the concept of localized magnons was used for a detailed analysis of the low-temperature magnetothermodynamics 
in the vicinity of the saturation field for two representative systems, the sawtooth chain (or A-chain) and the kagome 
lattice. In particular, the authors of these papers mapped the low-energy degrees of freedom of the sawtooth chain 
(the kagome lattice) to the hard-dimer gas on a one-dimensional lattice (the hard-hexagon gas on a triangular 
lattice) and used the results for the classical lattice gases to discuss the properties of the spin systems. They also 



provided exact diagonalization data for finite sawtooth chains (up to A^ = 20 sites) to illustrate the efficiency of 
the hard-dimer description of the spin chain at low temperatures near the saturation field. 

In the present paper we extend the previous studies on the low-temperature strong-field magnetothermody- 
namics examining various other frustrated quantum spin antiferromagnets supporting localized magnon states. We 
emphasize that all such spin systems exhibit a universal behavior. It is determined by a specific hard-core object 
lattice gas which mimics the independent localized magnon states. Thus, the one-dimensional hard-dimer behavior 
is also inherent in the two-leg ladder or kagome-like chains, whereas the hard-hexagon behavior is also inherent 
in the star lattice (see Table ^ below) . Moreover, we examine the thermodynamic properties of the spin models 
which are described by a gas of monomers. We also provide exact diagonalization data for finite spin systems 
of up to A^ = 30 sites and discuss to what extent a hard-core object description can reproduce the properties of 
the spin systems. We compare analytical and numerical results for different one-dimensional lattices and for the 
two-dimensional square-kagome lattice. 

The remainder of the paper is organized as follows. In Sec.|21we introduce the spin models to be discussed. In 
Sec. 121 we discuss a hard-core object description of the independent localized magnon states. We also consider how 
the lattice gas description can be extended for a wider temperature/field region. In Sec. 0] we present the exact 
diagonalization data for finite systems and compare them with the analytical predictions which follow from the 
hard-core object picture. We discuss briefly the effects of the localized magnons on the temperature dependence 
of the specific heat and the magnetocaloric effect near the saturation field. Finally, in Sec. |31 we summarize our 
findings. 

2 Geometrically frustrated spin models and localized magnons 

In this paper we consider several frustrated quantum spin lattices discussed so far in the literature by various 
authors, namely, the diamond chain [14], the dimer-plaquette chain [15], the sawtooth chain [16], the two-leg 
ladder [17], two kagome-like chains [18,19] (Fig. ^|, the square-kagome lattice [20], the kagome lattice [3,5,21], 
the star lattice [3,22], the checkerboard lattice [23] (Fig. |2J). (Note that further models hosting localized magnons 
can be constructed.) The Hamiltonian of the system consisting of N antiferromagnetically interacting spins in a 
magnetic field h reads 

H = y] J,,s, . s, -hS^^Y] J,, (I {s+s- + S-S+) + s^s^ - hS\ S^=y^s! (1) 



(y) (»j) ' 

with sf = s{s + 1) and Jij > 0. From Refs. [3-7] we know that for certain relations between Jy (see Tabled]) all 
spin lattices shown in Figs. ^ and |21 support localized magnon states which become relevant at strong magnetic 
fields. In what follows we fix the energy scale by setting either J = 1 (for systems with only one exchange integral) 
or Ji = 1 (for systems with more then one exchange integral). For the diamond and dimer-plaquette chains and 
for the two- leg ladder the magnons may be localized on the vertical bonds (see panels a, b and d in Fig. ^, for the 
sawtooth chain the localized magnons may be trapped in a valley between two neighboring triangles (panel c in 
Fig. ^, whereas for other lattices the localized magnons may occupy the even polygons shown by fat lines in Figs. 
^ 121 The explicit expressions for the localized one- magnon state having the smallest possible region of localization 
read 

|llm) ex (|s- li,S2) - |si,s- I2)) |s, . ..,s) 

(panels a, b and d in Fig. ^, 

|llm) ex (|s - li, S2, S3) - 2|si, s - I2, S3) + jsi, S2, s - I3)) |s, . . . , s) 
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Figure 1: (Color online) Several one-dimensional frustrated quantum spin lattices supporting localized magnon 
states: the diamond chain [14] (a), the dimer-plaquette chain [15] (b), the sawtooth chain [16] (c), the two-leg 
ladder (the spins are sitting only on the squares, not on the intersections of the diagonals) [17] (d), and two 
kagome-like chains [18,19] (e and f). The trapping cells occupied by localized magnons are shown by fat lines. 
Below each lattice we show the corresponding auxiliary lattice filled by hard-core objects (monomers (a and b) and 
dimers (c, d, e and f)). 
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Figure 2: (Color online) Several two-dimensional frustrated quantum spin lattices supporting localized magnon 
states: the square-kagome lattice [20] (a), the kagome lattice [3,5,21] (b), the star lattice [3,22] (c), and the 
checkerboard lattice (the spins are sitting only on the squares, not on the intersections of the diagonals) [23] 
(d). The trapping cells occupied by localized magnons are shown by fat lines. Below each lattice we show the 
corresponding auxiliary lattice (square (a and d) and triangular (b and c) lattices) filled by hard-core objects 
(monomers (a), hexagons (b and c) and squares (d)). 
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Table 1: Specific and generic properties of the spin systems shown in Figs. n]and[21 ei is the difference between 
the energy of the fully polarized state and the energy of the one-magnon state in zero field and determines the 
saturation field hi — ei, ri,„ax is the maximum number of independent localized magnons for the closest possible 
packing, M is the number of sites in the auxiliary lattice with hard-core objects, N is the number of sites in the 
spin lattice. 
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(the sum over m runs, say, counterclockwise over the sites ?Ti of a trapping polygon in panels e and f in Fig. Qand 
in all panels in Fig. ^. Here s^ or s — 1^ denote the values of sf and s is the maximal possible value of sf . We start 
with the zero-field case. The energy of the localized one-magnon state can be written as -Bfm — ei, where -Bfm is 
the energy of the fully (ferromagnetically) polarized spin system and the values of ei for the considered systems 
are given in Table Q] We can fill the spin lattice by n = 1, . . . , ?T-max, '^max oc N (see Table P) independent localized 
magnons. Each localized magnon decreases the magnetization S^ by 1 and the state with n localized magnons has 
S^ = Ns — n. The independent localized magnon states are the lowest-energy states in the corresponding sectors 
of S^ = Ns — 1, . . . , Ns — Timax [4, 24]. Although the localized magnon states exist also for anisotropic XXZ 
Heisenberg exchange interaction in (^ and arbitrary spin length s, in what follows we restrict ourselves to the 
isotropic (i.e. XXX) Heisenberg exchange interaction and (in numerical computations) to the extreme quantum 
case s — 1/2 without loss of generality for discussion of the universal low-temperature behavior of the spin systems 
near the saturation field. However, we have to bear in mind that the effects of the localized magnon states are true 
quantum effects which become less and less pronounced as s —^ oo [5, 11]. 

Consider now the localized magnon states in the presence of an external field h. A localized one-magnon state 
has the energy Ei{h) — E-pM ~ hsN — ei + h, whereas a state with n independent localized magnons has the energy 



En{h) = E'fm - hsN -n{ei-h). 



(2) 



Evidently, n independent localized magnons can be placed on a lattice with N sites in many ways. We denote by 
gN{n) > 1 the number of ways in which n independent localized magnons, each occupying the smallest possible 
area, can be put on a lattice with N sites. The energy of any of these gN{n) states is the same and is given by 
Eq. J^Jl, i.e. gN{n) is the degeneracy of the independent localized n-magnon states. At the saturation field hi ~ ei 



the energy En{hi) (j^J does not depend on n that further increases the degeneracy of the independent locahzed 
magnon states. 

Let us briefly recall the consequences of the localized magnons [3-13]. Due to the locahzed magnons different 
spin systems have identical/universal behavior in the ground state in strong magnetic fields. Since the lowest- 
energy states with different numbers of independent localized magnons n — 1, . . . , Timax have the same energy at 
the saturation field hi, E-pM ~ hisN, the ground-state magnetization exhibits a jump at h = hi between the values 
Ns and Ns — rimax- This jump is accompanied by a preceding plateau. Next, the ground-state energy at the 
saturation field exhibits a huge degeneracy which grows exponentially with the system size (at least as 2"°""=) and 
hence the ground-state entropy per site remains finite at h = hi. Finally, a lattice distortion which preserves the 
symmetry of the cell which hosts a localized magnon may lower the total ground-state energy (which consists of 
the magnetic and elastic parts). This deformation, which obviously cannot exist for h> hi, leads to a field-tuned 
ground-state structural instability in the vicinity of the saturation field. 

One may expect that the contribution of localized magnons is dominant also at low temperatures and strong 
magnetic fields around the saturation field. The analysis for the sawtooth chain and the kagome lattice reported in 
Refs. [8,12,13] supports this expectation. We now move on to discuss the contribution of the independent localized 
magnons to the low-temperature properties of the considered spin systems in the vicinity of the saturation field. 

3 Independent localized magnons and classical lattice gases 

In this Section, we examine the contribution of the independent localized magnons to the thermodynamic quantities 
of the spin systems. Moreover, we use for this purpose a hard-core object representation of the independent localized 
magnons of the smallest area. Such a representation allows us to utilize the broad knowledge on classical lattice 
gases. In the context of this approach a number of questions naturally arise which we have to consider. First, one 
has to show that the independent localized magnon states of smallest area are linearly independent in each sectors 
of S"^ — Ns — 1, ... , Ns — riniax- Only in this case they all contribute to the partition function of the spin system. 
Although there is no proof of the linear independence of the independent localized magnons in the general case, 
this question has been discussed in Refs. [8, 13] and a rigorous proof for all spin systems shown in Figs. miHis given 
in Ref. [25]. Thus, it was shown that for the diamond chain, the dimer-plaquette chain, the two-leg ladder, and the 
square-kagome lattice (the orthogonal type of frustrated spin systems in the nomenclature of Ref. [25]) as well as 
for the sawtooth chain and two kagome-like chains (the isolated type of frustrated spin systems) the independent 
localized magnon states are linearly independent in every sector of S^ — Ns — 1, . . . , Ns — rimax- Moreover, for 
a system of the orthogonal type they form an orthogonal basis. The kagome lattice, the star lattice and the 
checkerboard lattice belong to the frustrated spin systems of codimension one type. They possess exactly one 
nontrivial linear relation between their localized one-magnon states, but the set of n independent localized magnon 
states is linearly independent for n — 2, . . . ,ninax [25]. Thus, in the thermodynamic limit the linear independence 
of the localized magnons is fulfilled for all considered spin systems. 

Next, we have to check that there are no other low-energy states (except the independent localized magnons) 
in every sector of S^ — Ns — 1, ... , Ns — rimax, or at least to show, if they do exist, that their contribution is 
vanishingly small in the thermodynamic limit. In Refs. [8, 13] it was shown that for the kagome lattice there are 
indeed a few classes of states in the sectors S^ — Ns ~ 1, . . . , Ns ~ nmax which cannot be represented in terms 
of the independent localized magnons of smallest area. Thus, it remains unclear whether these additional states 
contribute to the partition function in the thermodynamic limit. In our numerical data (see Tables |21 and |3| below) 
we also find extra spin states for the diamond chain with J2 = 2Ji, the square-kagome lattice with J2 = Ji, the 
two- leg ladder with J2 = 2Ji and two kagome-like chains (but not for the diamond chain with J2 = 3Ji, the 
dimer-plaquette chain with J2 = Ji, J3 = 2Ji, the square-kagome lattice with J2 = 2Ji, the sawtooth chain and 



the two- leg ladder with J2 = 3Ji). For all these systems, except the square-kagome lattice with J2 = Ji, the 
numbers of the extra spin states do not increase with increasing of the lattice sizes that indicates their irrelevance 
in the thermodynamic limit. For the square-kagome lattice with J2 — Ji, as for the checkerboard [25] and the 
kagome [26] lattices, we find much more extra spin states, however, our data restricted to finite systems are not 
sufficient to find reasonable tendencies for the thermodynamic limit. Furthermore, we have to discuss whether the 
independent localized magnon states are separated by a finite energy from higher-energy states. We postpone a 
discussion of this issue till Sec. ^ where the results of the exact diagonalization for finite systems are presented. 

We are interested in the contribution of the independent localized magnons of smallest area to the thermody- 
namics of the spin systems. To find this contribution we must calculate the part of the partition function of the 
spin system, Z{T, h, N) — ^ exp {~Ej{h, N)/kT) {j denotes all states of the system of N spins s on a lattice), 
which comes from the independent localized magnons. Denoting this quantity by Z\y^{T, h,N) we have 

ZUT,h,N) = '5; 5.(n)exp (~M^) = exp (~^^^^r^) E ff^^exp {^n) (3) 

n=0 \ / \ / „^o 

with fi ~ ei — h. 

An important step for further calculation of Zini(r, h, N) (jS)) is a mapping onto a hard-core object lattice gas. 
Consider an auxiliary lattice in which each site can be occupied or not by a hard-core object (monomer, dimer, 
hexagon or square) which corresponds to a localized magnon of smallest area (see Figs. n]|21l. Let us denote by A/" 
the number of sites of the auxiliary lattice; the relation between N and M for the considered spin systems is given in 
Tabled We note that gN(n) is simply the canonical partition function Z{n,J\f) of n hard-core objects placed on a 
lattice with J\f sites. Furthermore, 'E.{T, fi,J\f) = '^^exp{fj.n/ kT) Z (n, J\f) is the grand canonical partition function 
of hard-core objects placed on a lattice with A/" sites and fj, is the chemical potential. As a result, according to Eq. 
lO we arrive at the basic relation between the thermodynamic quantities of independent localized magnons and 
hard-core models, i.e. Zini(T, h, N) = exp(-(£;FM - hsN)/kT)E{T, ^,7V). 

The calculation of S(T, /x. A/") for a classical hard-core object lattice gas (usually in the thermodynamic limit), 
in general, is a nontrivial problem [27]. ft may be convenient to perforin such calculations using the occupation 
number representation. Let us introduce the occupation number rii which takes two values, and f, depending 
whether the site i of the auxiliary lattice is empty or occupied, i = 1, . . . ,J\f. Then the grand canonical partition 
function of hard-core objects can be rewritten in the following way 

E{T,^l,^f)= E •■• E exp(^EnJi?(ni,...,nA,), (4) 

ni=0,l nj^=0,l \ 1=1 / 

where the function R{ni, . . . ,nj\f) ^ I appears if the hard-core objects are extended (dimers, hexagons, squares) 
rather than simple monomers. For example, for the sawtooth chain, for which the corresponding hard-core objects 
are rigid dimers, we have R{ni, . . . , n_\f) = J|,-(f — riirii+i) and all the terms in Q which correspond to two adjacent 
sites of the auxiliary lattice being occupied have dropped out. 

Using the relation between the independent localized magnons partition function Zi„^{T,h, N) and the grand 
canonical partition fmiction of hard-core object lattice gas S(T, /i, A/") we get the following result for the Helmholtz 
free energy (per site) of the independent localized magnons 

Fin-,iT,h,N) ^E^M_j^^_ ^^^ lnS(r,/i,AA) 

The entropy Si^{T,h,N) = ~dFy^{T,h,N)/dT and the specific heat Ci^{T,h,N) = TdSu^{T,h, N)/dT follow 
immediately from Eq. (0. We can also calculate the average number of hard-core objects (n) = J2i=ii''^i) ~ 
kTdlnE{T, ^,Af)/diJ, which yields the magnetization {3")^^ — sN — (n). 

Next, we turn to the specific lattice gases which are related to the considered spin systems. Although the results 
for various lattice gases are not new and can be found in the literature [27] we present some of them here for easy 
references in view of further discussions for the spin systems. 
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3.1 Monomers 

The independent localized magnon states in the diamond and dimer-plaquette chains as well as in the square- 
kagome lattice can be mapped onto a lattice gas of monomers (see panels a and b in Fig. ^ panel a in Fig. |21 
and Table ^. The calculation of the thermodynamic quantities for the lattice gas of monomers is simple. We have 
Z{n,N) - CJ^ = m/{n\{N -n)\), E{T,fi,Af) = T,n=o ^J^/ ^Mf^n / kT) = (exp(/i//cT) + 1)-^ or, using Eq. gj, 
E{T,iJ,J\f) = {J2n=o 1 exp(/in/fcr))-'^ = (1 + cxp{ii/kT))-^ . The Helmholtz free energy of the gas of monomers 
reads 

Fin.{T,h,N) EpM , AA / ^i^ 



and therefore 



= T7 li^l + expx) - — — , (7) 



kN N \ 1 + exp X 

CUT,h,N) Af (f) 



2 

(8) 



kN N cosh^ f ' 

{S^)im _ lA/" expcc 
sN ~ s N 1 + expx ^ ^ 

with X — ^/kT. Eqs. lO, ^ are symmetric with respect to a; ^ —x. We note that at saturation /i = and Eq. 
iJTIl yields S'im(T, /ii, N)/kN = (Af/N) In 2. That is the value of the residual entropy per site obtained earlier for 
the diamond and dimer-plaquette chains [11]. The specific heat Ci^{T, h,N)/kAf, Eq. ||SJ), exhibits two symmetric 
maxima of height C„iax ~ 0.43922884 at a; « ±2.39935728. The universal dependences of (iVs- (5^)im)/7V, Eq. Q, 
of S'lm/fcA/', Eq. 0, and of Cim/fcA/", Eq. ((SJ, are shown in Figs. |31 El (for a more detailed discussion of these Figures, 
see Sec. EI). 

We may consider an improvement of a lattice gas description of low-lying energy levels of spin systems taking 
into account higher-energy states which are separated from the hard-core object states by a gap. According to Eq. 
Q within a lattice gas model there are only two states at each lattice site. Namely, the site may be either empty 
(rii = 0) and such a state occurs with the probability ex 1, or the site may be occupied by monomer (n^ = 1) and such 
a state occurs with the probability oc exp(/i//cT). We may increase the number of states at each sites introducing 
an additional higher-energy state of monomers which occurs with the probability oc exp((/i— U)/kT) and U is large 
positive on-site energy. Obviously, if [/ -^ oo the higher-energy states are irrelevant, but for finite U they become 
relevant as temperature grows. In the spin picture, the second state corresponds to a higher-energy state with a 
wave function located in the area which is mapped onto one site of the auxiliary lattice. Our calculations confirm 
that the made assumption really yields an improvement for higher temperatures (see the corresponding panels in 
Figs.OllZllHI). The grand partition function of the modified lattice-gas model can be easily calculated 

S(r,/i,AA)=(^l + exp^+exp^) . (10) 

From Eq. H10|) we find the entropy, the specific heat, and the magnetization using usual thermodynamic relations. 
The corresponding results are illustrated in Figs. |2| 13 El assuming values for U related to the energy separation 
between the localized magnon states and higher-energy states for finite spin lattices, see Sec. ^ The improved 
approximation agrees with exact diagonalization results up to higher temperatures in comparison with the monomer 
lattice gas predictions © - 113 . 
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Figure 3: The universal dependences of {sN — {S^)\m)/J^ ©j Sim/kAf and C\ni/kJ\f © (monomer universality 
class) together with the corresponding results for finite spin systems (diamond chain with A^ = 18, Ji = 1, J2 = 3 
[left column], dimer-plaquette chain with N = 20, Ji = J2 = 1, J3 = 2 [middle column], square-kagome lattice 
with N ~ 18, Ji = 1, J2 = 2 [right column]). The universal dependences are shown by solid lines, the exact 
diagonalization data are shown by symbols. Note that the exact diagonalization data for the lowest temperature 
coincide with the universal dependences for monomers. We also show by thin broken lines the analytical results 
which follow from pU|l with U = 1.57, 0.58, 1.34 for the diamond chain, dimer-plaquette chain and square-kagome 
lattice, respectively, which are in good agreement with numerical data for higher temperatures. 
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Figure 4: The universal dependences for monomers (solid lines) of {sN— (5^)im)/A/', S\m/kN and C\m/kM together 
with the corresponding results for finite spin systems belonging to the monomer universality class but which have 
extra (not hard-monomer) spin states: the diamond chain with TV = 18, Ji = 1, J2 = 2 (left column), the square- 
kagome lattice with N — 18, Ji = J2 = 1 (right column). The numerical results for spin systems are shown for 
temperatures kT = 0.01 (filled symbols) and kT = 0.1 (empty symbols). 
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3.2 One-dimensional dimers 

The independent localized magnons in the sawtooth chain, the two-leg ladder and two kagome-like chains can be 
mapped onto a one-dimensional lattice gas of rigid dimers (see panels c, d, e, f in Fig. ^and Table ^. The grand 
partition function S(T, /ijA/") Q of rigid dimers on a chain of Af sites with periodic boundary conditions can be 
calculated using the transfer-matrix method. The resulting expression for S(T, /ijA/") reads 



\Ar , \Ar \ 1 , /I , ._.„ ^^ 



S(T,/i,AA) = A^+A^, Ai,2 = -±^--l-exp^. (11) 

In the thermodynamic limit Af —^ oo only the larger eigenvalue of the transfer matrix plays role in l|ll(l and we 
have 



FiUT,h,N) Efm , A/".^, /I /I ^i \ ,.^. 

For the entropy, specific heat, and magnetization we find from Eq. H12() 




Mz^4,„^w^,_i_.i, ___,,, 



4\/i+expa:;, 



Ci„,{T,h,N) TV x^expa; 



kN iVg^i^^^p^) 



iS' 




'^ -^^^2 , /TT^,^ 



(14) 



(15) 



respectively. Note that Eqs. ifT^ . !fT^ are not symmetric with respect to x ^ —x. Eqs. ((T^ . ffT^i . ((T3|l 
were examined in Refs. [8, 12, 13] in the context of the sawtooth chain. We also note that at saturation, i.e. 
at /i = 0, Eq. H13|l gives the values of the residual entropy per site for the kagome-like chains and the two- leg 
ladder, Sin,{T,hi,N)/kN = {M/N)\n{{l + V5)/2) [11]. The universal dependences {Ns - (5^)inO/A/', Eq. |T^. 
S\m/kM, Eq. p3|l . and Cim/kAf, Eq. p4|) . are shown in Figs. [S] and (for a more detailed discussion of these 
Figures, see Sec. lU. The specific heat Cim(T, h, N)/kM exhibits two maxima of heights C,'^|',J* « 0.34394234 and 
&^^^ « 0.26887020 at [h - hi)/kT « 2.81588498 and {h - hi)/kT « -4.05258891, respectively 

Again we may consider an improvement of a lattice gas description of the low-lying energy levels of spin 
systems. One way to consider higher-energy states was suggested by M. E. Zhitomirsky and H. Tsunetsugu [8]. 
These authors treat the factor 1 — riin^+i in Q), as the F ^- oo limit of a generalized factor eyi^{— {V / kT)nini+i) = 
1 + {exp{—V/kT) — l)nini_|_i. We may relax the hard-core restriction assuming the intersite interaction V to be 
large but finite and approximately equal to the energy separation between the localized magnon states and higher- 
energy states (in Ref. [8] the energy gap was also estimated using variational calculations). As a result, we arrive 
at a one-dimensional lattice gas with a finite nearest-neighbor repulsion. The grand partition function of such a 
gas is calculated by the transfer-matrix method and reads 

A/" I \A^ 



E{T,ti,Af) - Af+A: 



2 ? 



11 fi-V I n 1 /i-y 1 2{p-V) 

Ai 2 = - + - exp , ^ ± \ - + exp -^ exp , ^ + - exp — ^ -. (16) 

2 2 ^ kT \ A ^ kT 2 ^ kT A ^ kT ^ ' 

The thermodynamic quantities which follow from Eq. (|16|l are in a good agreement with exact diagonalization 
data for the sawtooth chain up to higher temperatures in comparison with the predictions in the V^ — > cx) limit [8] . 
In the spin picture a finite V corresponds to higher-energy states with wave functions located on the areas which 
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Figure 5: The universal dependences of {sN— (5^)im)/A/', Sim/kM and Cim/kAf (one-dimensional dimer universality 
class) together with the corresponding results for finite spin systems (sawtooth chain with N — 20, Ji — 1, J2 = 2 
[left column] and two-leg ladder with N = 20, J = 1, J2 = 3 [right column]). The universal dependences which 
follow from Hll|) are shown by thick solid {JV -^ 00) and dashed lines {JV = 10), the exact diagonalization data are 
shown by symbols. Note that the exact diagonalization data for the lowest temperature coincide with the universal 
dependences for dimers. We also show by thin broken lines the analytical results which follow from (|18() for A/" = 10 
with U = 1.00 and 1.10 for the sawtooth chain and two-leg ladder, respectively, which are in good agreement with 
numerical data for higher temperatures. 
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Figure 6: The universal dependences of {sN — (5^)im)/A/', S\m/kM and C\ni/kJ\f for onc-dimcnsional dimers 
(A/" ^ cx) - thick solid lines; finite Af - dashed lines) together with the corresponding results for finite spin systems 
belonging to the one-dimensional dimer universality class but which have extra (not hard-dimer) spin states: the 
two-leg ladder with A^ = 20 (corresponds to JV = 10), Ji = 1, J2 = 2 (left column), the kagome-like chain shown 
in Fig. ^ with A^ = 18 (corresponds to M — 6), J — 1 (middle column) and the kagome-like chain shown in Fig. 
^ with A^ = 20 (corresponds to A/" = 4) , Ji = 1 (right column) . The numerical results for spin systems are shown 
for temperatures kT = 0.01 (filled symbols) and kT = 0.1 (empty symbols). 
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are mapped onto two adjacent sites of the auxiliary one-dimensional lattice. In addition, we can consider a similar 
improvement as for the monomer problem, namely we can assume additional higher-energy states of dimers which 
occur with the probability oc exp((/x — U)/kT), where U is positive and large. That means that more than one 
dimer is allowed on the same site, but that cost additional energy U. In the spin picture, we take into account 
a higher-energy state with a wave function located in the area which is mapped onto one site of the auxiliary 
one-dimensional lattice. Thus, combining the assumption about higher-energy state of dimer at each site and the 
weakening of the hard-core restriction we arrive at the following formula for the grand partition function 



E{T,h,M) = Xf + X^ + X: 



3 ' 



H-V\2 „..„-" /^ -,,„..„ /^ U\ _^f yW^. 2^i-U-V 



K' 


+ Af, 








I- 


^i + exp 


kT 


+ exp 


kT ' 



A^-^l + cxp:i^jA^-exp^^l + exp^--j-exp^--jjA + exp^^3^^.0; (17) 

in the thermodynamic limit A/" — » cxd only the largest root of the qubic equation enters the thermodynamic quan- 
tities. In the limit U ^ oo Eq. H17|l transforms into Eq. (|16|l . In what follows we consider only the other limit 
y -^ oo when Eq. (|17|l becomes 

EiT,^l,^f) 

From Eq. (|18|l we find in a usual way the entropy, the specific heat, and the magnetization, which are shown in 
Figs. IHl assuming values for U again obtained by inspection of the energy spectrum of finite spin lattices (for 
more details see Sec.EJ. 

3.3 Hexagons on a triangular lattice and squares on a square lattice 

For completeness we give here some remarks on the two-dimensional lattice gases of hard-core objects which emerge 
in the context of localized magnons in some two-dimensional frustrated quantum antiferromagnets. Such lattice 
gases are much more complicated subject for rigorous analytical or numerical studies. 

For the kagome and star lattices the independent localized magnons of smallest area can be put in correspondence 
to the hard hexagons on a triangular lattice (see panels b and c in Fig. [21and Table^. The grand canonical partition 
function for the hard hexagons on triangular lattice, S(T, /i. A/"), was found by R. J. Baxter (hard- hexagon problem, 
see Ref. [27]). The parametric dependence of (S(T, /x, A/"))^/ on the activity z — exp{fi/kT) in the thermodynamic 
limit J\f -^ oo can be found in Ref. [27]. The peculiarity of the hard-hexagon model is the existence of a phase 
transition at a critical value of activity z^ = (U + 5V5)/2 w 11.09016994 {xc = In^c « 2.40605913). The low- 
activity phase is uniform (liquid) with identical average occupation numbers of hexagons on each of three triangular 
sublattices, ni = n2 = n^, whereas the high-activity phase is nonuniform (solid) with rii > n2 — n^. The value 
of the entropy per site at z = 1 is 0.33324272fc; this number gives the residual entropy per site due to localized 
magnon states of smallest area for the kagome and star lattices [8, 11, 13] . The hard-hexagon problem in the context 
of the low-temperature magnetothermodynamics of the kagome lattice in the vicinity of the saturation field was 
discussed in detail in Refs. [8, 13]. 

For the checkerboard lattice the independent localized magnons of smallest area can be put in correspondence to 
the hard-core squares on a square lattice, where the size of the hard squares corresponds to the nearest-neighbor and 
the next-nearest-neighbor exclusions (see panel d in Fig. Inland Tabled]). We are not aware for an exact analytical 
solution of such a model, but approximative results are available [28]. In the context of the low-temperature 
magnetothermodynamics of the checkerboard lattice in the vicinity of the saturation field the hard-square problem 
has been discussed recently in Ref. [29] . The value of the entropy per site dX z — 1 obtained applying a classical 
Monte Carlo method is about 0.2946fc [29] (a direct calculation for periodic 8 x 8 (8 x 10) lattice yields about 
0.2949fc (0.2948fc) [30]). 
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4 Exact diagonalization versus hard-core object description 

After presenting explicit expressions for the thermodynamic quantities by using the correspondence between the 
locahzed magnon states of the quantum spin system and a classical hard-core description we now test our analytical 
results by comparison with finite-lattice numerical results of the full spin-1/2 system. Note that the classical 
equations for the monomer problem, see Sec. 13. II do not depend on the size of the system, whereas the corresponding 
classical equations for the dimer problem, see Sec. 13.21 are size dependent. For the number of degenerate states 
presented in Tables |21 and 13 we could consider systems of up to A^ = 30 spins, since we can restrict the calculation 
of sectors with high total S^. For the thermodynamics we used the data of full diagonalization of spin systems 
of either A^ = 18 or A^ = 20. For the estimation of the temperature and field ranges, in which the hard-core 
description is valid, and also for the estimation of the on-site energy parameter U, see Secs. l3.ll and l3.2l it is useful 
to find a measure for the thermodynamically relevant energy separation A between the localized magnon states 
and the other eigenstates of the spin system. For that we have calculated the integrated low-energy density of 
states at saturation field and define A as that energy value above the localized magnon ground-state energy, where 
the contribution of the higher-energy states to the integrated density of states becomes as large as the contribution 
of the locahzed-magnon states. 

In Figs. 121 and 01 we present the results for the magnetization {Ns — {S^))/Af, the entropy S/kAf, the specific 
heat C/kAf in dependence on {h — hi)/kT for spin systems belonging to the monomer universality class, i.e. for 
the diamond chain, the dimer-plaquette chain and the square-kagome lattice. Note that the universal formulas 
© - @ depend only on {h — hi)/kT, i.e. they are identical for all temperatures kT — 0.01, 0.1, 0.2, 0.3, 0.5 
considered in Figs. and ^ 

Fig. corresponds to spin systems, for which the number of degenerate spin states at the saturation Wsp equals 
the number of states of the corresponding hard-core object lattice gas W, see Table (21 Furthermore, these spin 
systems exhibit a quite large energy separation A. We find A « 1.5 for the diamond chain (Ji = 1, J2 — 3, 
A^ = 18), A « 0.6 for the dimer-plaquette chain (Ji = J2 = 1, J3 = 2, A^ = 20), and A « 1.5 for the square-kagome 
lattice (Ji = 1, J2 = 2, A^ = 18). Obviously, in Fig. the curves for the spin systems at temperatures up to 
kT « 0.1 « O.lA are almost identical with the universal curves. But also for kT > 0.1 the qualitative behavior 
of the entropy, the specific heat, and the magnetization is quite well modeled by the universal formulas Q, Q, 
0. The modified lattice-gas model, see Eq. (|10|l . with parameters U related to A indeed leads to an improved 
quantitative agreement with the results for the spin systems at higher temperatures. 

Fig. 21 corresponds to spin systems, for which the number of degenerate spin states at the saturation field Wsp is 
larger than the number of states of the corresponding hard-core object lattice gas W, see Table|21 However, for the 
diamond chain with J2 = 2 Ji, we have only one extra spin state, Wsp = W + 1, which becomes irrelevant for larger 
A^. By contrast, there is a noticeable disagreement between Wsp and W for the square-kagome lattice with J2 — Ji 
(Table I2I1. Similarly to the kagome [8,13,26] and the checkerboard [25] lattices, for the square-kagome lattice 
with J2 — J I we have extra spin states which are not covered by the hard-core description illustrated in Sec. 13. l| 
The number of these extra spin states depends on system size A^. From our numerical data for A^ = 18, 24, 30 
we are not able to conclude whether the localized magnon states of smallest area dominate the low-temperature 
thermodynamics in the thermodynamic limit. Moreover, for these systems with extra spin states we estimate 
significantly lower separations A of higher-energy states, namely A « 0.3 for the diamond chain (Ji = 1, J2 = 2, 
A^ = 18) and A w 0.2 for the square-kagome lattice (Ji = 1, J2 = 1, A = 18). As a consequence we may find (see 
Fig. ^ a perfect agreement between the spin and the hard-core data only for very low temperatures, but realize a 
quantitative deviation already for kT w 0.1. For the square-kagome lattice with Ji = 1, J2 = 1 the larger number 
of extra spin states leads to a more pronounced deviation from the universal formulas. Nevertheless, the numerical 
data illustrate that the localized magnon states of smallest area contribute substantially to the partition function 
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Table 2: Number of degenerate states: lattice gas of monomers and spin systems. For a gas of monomers on A/" 
sites we have W = S(T, ^ = 0,7V) = 2^. 



M 


W 


diamond 

(J2 = 2Ji) 


diamond 

(J2 = 3Ji) 


dimer-plaquette 
{Js = 2Ji) 


square-kagome 

(J2 = Jl) 


square-kagome 
(J2 = 2Ji) 


3 


8 


9 


8 


8 


18 


8 


4 


16 


17 


16 


16 


42 


16 


5 


32 


33 


32 


32 


52 


32 


6 


64 


65 


64 


64 


- 


- 


7 


128 


129 


128 


- 


- 


- 


8 


256 


257 


256 


- 


- 


- 



Table 3: Number of degenerate states: one-dimensional hard-dimer lattice gas and spin systems. For one- 
dimensional dimer gas on ^f sites we have W = E{T, /i = 0,7V) = Ai(^ = 0)-^ + A2(/i = 0)^. Since Ai(^ = 0) = $ 
is the golden number and X2{fJ- = 0) = 1 — $ = —1/$ we may use the properties of powers of $ and its recipro- 
cal. Namely, for any even integer n <&" -f- 1/$" is a whole number; in particular, 2, 3, 7, 18, 47, 123, 322 for 
n = 0, 2, 4, 6, 8, 10, 12. (For any odd integern $"-l/i>" is a whole number; in particular, 1, 4, 11, 29, 76, 199 
for n = 1, 3, 5, 7, 9, 11.) See: http://goldennumber.net/phipower.htm 



A^ 


W 


sawtooth 


two-leg ladder 
(J2 = 2Ji) 


two-leg ladder 
(J2 = 3Ji) 


kagome-like. 
Fig. IB 


kagome-like, 
Fig.Cf 


4 


7 


7 


8 


7 


9 


9 


6 


18 


18 


19 


18 


20 


20 


8 


47 


47 


48 


47 


49 


- 


10 


123 


123 


124 


123 


- 


- 


12 


322 


322 


323 


322 


- 


- 



thus leading to a low-temperature behavior which is at least qualitatively in agreement with the universal behavior 
given by Eqs. Q, Q, and I©. 

Next, we discuss the results for the magnetization (TVs— (5^))/7V, the entropy S/kAf and the specific heat C/kJ\f 
in dependence on {h — hi)/kT for spin systems belonging to the dimer universality class, i.e. for the sawtooth chain, 
the frustrated two-leg ladder and the kagome-like chains of type I and II, see Figs. [S] and Again, the universal 
formulas Hll|l - (|15|l depend only on {h—hi)/kT, i.e. they are identical for all temperatures kT = 0.01, 0.1, 0.2, 0.3 
considered in Figs. |31 and IHl However, as mentioned already above, they are size dependent. Hence, we compare the 
spin data for systems of finite size TV with the hard-core data of corresponding size 7V (cf. Table ^ obtained from 
Eq. Hll|l. Note, however, that the systems considered in Fig. [S] and in Fig. (HJ left panel, correspond to 7V = 10, 
where the curves for the finite hard-dimer system are already close to those for the thermodynamic limit. 

The results shown in Fig. [S] belong to systems with no extra spin states, i.e. Wsp — W, see Table 13 and quite 
large energy separation A w 0.8 (sawtooth chain, N — 18) and A w 0.9 (two-leg ladder, Ji == 1, J2 = 3, A^ = 20). 
Consequently, we find an excellent agreement between the spin data and the hard-dimer data for temperatures 
up to kT « 0.1 « 0.1 A. For larger temperatures again the qualitative behavior is well described by the universal 
formulas derived from Eq. Hll|l . The modified hard-core description based on Eq. (18) with appropriate parameters 
U leads to a further quantitative improvement, see Fig.|S| 
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Fig. El corresponds to spin systems, for which the number of degenerate spin states Wsp at the saturation field 
for finite spin systems is larger than the number of states W of corresponding finite hard-core object lattice gas, 
see Table 13 However, we have only one (two- leg ladder, Ji = 1, J2 = 2) or two (kagome-like chains I and II) extra 
spin states, which become irrelevant for larger N. The energy separation A for these systems with extra spin states 
we estimate to A « 0.9 for the ladder (Ji = 1, J2 = 2, TV = 20), A « 0.4 for the kagome-like chain I (N = 18), 
and A w 0.2 for the kagome-like chain II {N — 20). Therefore we obtain (see Fig. |SJ) a perfect agreement between 
the spin and the hard-core data only for very low temperatures, but realize a quantitative deviation already for 
kT « 0.1. Note that for the the kagome-like chain II, the size iV = 20 of the spin system corresponds to only 
A/" = 4, and the finite size effects are clearly largest in this case. 

From the experimental point of view the reported results manifest themselves most interestingly in a drastic 
change of the low-temperature specific heat, when the magnetic field passes the saturation field, and in the maximum 
of the isothermal entropy at saturation field leading to an enhanced magnetocaloric effect (for a general discussion 
of the magnetocaloric effect for quantum spin systems, see Ref. [12]). We illustrate that in Figs. [7| and |H1 where we 
present finite-lattice data for spin systems, but hard-core results for the thermodynamic limit to demonstrate that 
the discussed effects remain relevant for N —> 00. 

In Fig. \7\ we present the temperature dependence of the specific heat at three values of the external magnetic 
field, h = O.95/11, hi, l.Obhi, for several spin systems. Firstly, we note that within the hard-core object description 
the specific heat equals to zero at saturation field h = hi. That is also observed in a certain range of low 
temperatures around zero for the diamond chain with J2 > 2Ji, the square- kagome lattice with J2 > Ji or the 
two- leg ladder with J2 > 2Ji, but we do not observe the zero- value region of C/kN at low temperatures for the 
spin systems having extra spin states and low separation A. Secondly, we see that the specific heat in the vicinity 
of the saturation field (but h ^ hi) has an extra low-temperature peak which is satisfactorily reproduced within 
the hard-core object lattice gas approach (compare lines and empty triangles and circles in Fig. (T)). The hard-core 
description can be improved assuming a finite on-site energy parameter U (filled triangles and circles in Fig. [TJ. 
As it was already noticed in Ref. [8] in the context of the sawtooth chain and the kagome lattice, the experimental 
observation of the low-temperature peak of the specific heat can be a sign of highly degenerate localized magnon 
states. 

Finally, we consider in Fig. [Sladiabatic cooling processes, i.e. we show curves of constant entropy as a function of 
magnetic field and temperature. Again the exact diagonalization data for the diamond chain with N = 18, Ji = 1, 
J2 = 3 (left panel) and the two- leg ladder with N = 20, Ji = 1, J2 = 3 (right panel) are in a good agreement 
with the universal curves for monomers and one-dimensional dimers. Taking the data for spin systems, which show 
some asymmetry between h > hi and h < hi, we observe, that an efficient cooling to very low temperatures can 
be achieved by an adiabatic demagnetization via lowering the magnetic field from above saturation till saturation. 

5 Summary 

We have presented a universal description of thermodynamic properties of a wide class of frustrated quantum spin 
antiferromagnets at low temperatures in the vicinity of the saturation field. The reason for that is a dominant 
contribution of highly degenerate localized magnon states to the partition function, which can be calculated using 
a hard-core object representation. We have provided exact diagonalization data to justify such a picture. 

The spin systems hosting localized magnons can be grouped into different universality classes of hard-core lattice 
gases, namely gases of hard squares, hard hexagons, hard dimers or hard monomers. Comparing the analytical 
predictions with numerical data for finite spin systems we have found that the hard-core picture describes accurately 
the spin physics near the saturation field and at sufficiently low temperatures. Since other states of the spin 
system, not described by the hard-core model, become relevant as temperature grows, the hard-core description is 
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Figure 7: The temperature dependences of the specific heat C/kN vs. kT/hi for several finite spin systems. 

Monomer universaUty class: left panels, from bottom to top, the diamond chain with Ji = 1, J2 = 3, A^ = 18, the 

diamond chain with Ji = 1, J2 = 2, A^ = 18, the square-kagome lattice Ji — 1, J2 ~ 2, N = 18. One-dimensional 

dimer universality class: right panels, from bottom to top, the two-leg ladder with Ji = 1, J2 = 3, TV = 20, the 

two-leg ladder with Ji = I, J2 = 2, N = 20, the kagome-like chain shown in Fig. ^ with J = 1, A^ = 18. The 

thin solid lines correspond to h = 0.95/ii, the solid lines correspond to h = hi, the thick solid lines correspond to 

h = 1.05ft.. For comparison we present analytical results for monomers (left panels) and dimers. A/" — > 00 (right 

panels) (empty symbols; triangles correspond to h = 0.95fti, squares correspond to h ^ hi, circles correspond to 

h = 1.05ft.i). We also show by filled symbols the temperature dependence of the specific heat as it follows from 

Eq. (|10|) with U = 1.57 for the diamond chain and U = 1.34 for the square-kagome lattice or from Eq. (|18|l for 

Af —^ 00 with U = 1.10 for the two-leg ladder. 
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Figure 8: Constant entropy curves as a function of magnetic field and temperature, S{T,h,N)/kN = const, for 
the diamond chain with A^ = 18, Ji = 1, J2 = 3 (left) and the two-leg ladder with N = 20, Ji = 1, J2 = 3 (right). 
From bottom to top: S{T,h,N)/kN = 0.05, 0.10, 0.15, 0.20, 0.25. Solid lines correspond to analytical resuhs 
which follow from Eq. (Q (left panel) or from Eq. H13|) (right panel); broken lines correspond to analytical results 
for S{T, h, N)/kN = 0.20, 0.25 which follow from Eq. ^ with U = 1.57 (left panel) or from Eq. lP|) for TV ^ 00 
with [/ = 1.10 (right panel); symbols correspond to exact diagonalization data for finite spin systems. 

less accurate, but it remains qualitatively correct at higher temperatures. The hard-core object lattice gas models 
may be improved by relaxing the hard-core constraint. Such an improvement breaks the universality, but provides 
a better quantitative agreement with exact diagonalization results for higher temperatures and larger deviations 
from the saturation field. 

We emphasize that some peculiarities of the thermodynamic quantities arising due to the localized magnon 
states are of interest from the experimental point of view. Thus, the specific heat at the saturation field remains 
almost zero below a certain temperature whereas at fields slightly above/below the saturation field the specific heat 
exhibits a well-pronounced extra low-temperature peak. A frustrated spin system hosting localized magnons also 
exhibits a large magnetocaloric effect in the vicinity of the saturation field; similarly to an ideal paramagnet. 
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